Mon. Not. R. Astron. Soc. OOP. [TlfTil (20081 Printed 6 May 2008 (MN style file v2.2) 



The evolved stars of Leo II dSph galaxy from near-infrared 
UKIRT/WFCAM observations 

M. Gullieuszik^*, E. V. Held\ L. Rizzi^, L. Girardi\ P. Marigo^ and Y. Momany 

^ INAF/Osservatorio Astronomico di Padova, vicolo dell'Osservatorio 5, 1-35122 Padova, Italy 

^ Joint Astronomy Centre, 660 N. A'ohoku Place, University Park, Hilo, HI 96720, USA 

^ Dipartimento di Astronomia, Universita di Padova, vicolo dell'Osservatorio 2, 1-35122 Padova, Italy 

Accepted . . . Received . . . ; in original form . . . 



ABSTRACT 

We present a study of the evolved stellar populations in the dwarf spheroidal galaxy 
Leo II, based on JHKg observations obtained with the near-infrared array WFCAM 
at the UKIRT telescope. Combining the new data with optical data, we derived pho- 
tometric estimates of the distribution of global metallicity [M/H] of individual red 
giant stars from their V — Kg colours. Our results are consistent with the metallicities 
of RGB stars obtained from Ca ii triplet spectroscopy, once the age effects are consid- 
ered. The photometric metallicity distribution function has a peak at [M/H] = —1.74 
(uncorrected) or [M/H] = —1.64 ±0.06 (random) ±0.17 (systematic) after correction 
for the mean age of Leo II stars (9 Gyr). The distribution is similar to a Gaussian 
with (T[M/H] = 0.19 dex, corrected for instrumental errors. We used the new data to 
derive the properties of a nearly complete sample of asymptotic giant branch (AGB) 
stars in Leo 11. Using a near-infrared two-colour diagram, we were able to obtain a 
clean separation from Milky Way foreground stars and discriminate between carbon- 
and oxygen-rich AGB stars, which allowed to study their distribution in i^s-band 
luminosity and colour. We simulate the JHKg data with the trilegal population 
synthesis code together with the most updated thermally pulsing AGB models, and 
using the star formation histories derived from independent work based on deep HST 
photometry. After scaling the mass of Leo II models to the observed number of upper 
RGB stars, we find that present models predict too many 0-rich TP-AGB stars of 
higher luminosity due to a likely under-estimation of either their mass-loss rates at 
low metallicity, and/or their degree of obscuration by circumstellar dust. On the other 
hand, the TP-AGB models are able to reproduce the observed number and luminosi- 
ties of carbon stars satisfactorily well, indicating that in this galaxy the least massive 
stars that became carbon stars should have masses as low as ^ 1 Mq. 

Key words: Galaxies: individual: Leo II - Galaxies: stellar content - stars: AGB and 
post-AGB - stars: carbon - Local Group 



1 INTRODUCTION 

Leo II is one of the most distant dwarf spheroidal (dSph) 
satellites of the Milky Way. A number of pho tometric studies 
derive d quite different distances for Leo II. iMighell fc RichI 
1 19961 ). using the V magnitude of the horizontal branch 
(HB), placed the galaxy at a distance modulus (m — M)o = 
21.55 ± 0.18. U sing the I band magnitude of t he tip of the 
RGB fTRGBI. lBellazzini. Gennari fc Ferrarol (|2005i ) found 
(m — M)o = 21.84 ± 0.13. Distance estimates in the liter- 
ature are intermediate between these values, with typical 
uncertainties of ^ 0.2 magnitudes. 
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The estimates of mean metallicity o f Leo II range 
from [Fe /H] = -1-6 (IMighell fc RichI [19961 ) up to [Fe/H] 
~ —1.1 (|Dolphinl I2OO2I ). Recently, two independent spec- 
troscopic studies, based on the Call triplet method, 
found Leo II to be relatively met al-poor: a mean v alue of 
[Fe/H] = -1.74 was derived by iKoch et al.l (|2007l ) while 
iBosler. Smecker-Hane fc Stetson! (|2007l ) found [Fe/H] = 
-1.59. 

The Leo II dSph was considered for a long time as a 
typical "old" dSph. Then IMighell fc RichI (|l996l ) obtained 
an HST/WFPC2 colour magnitude diagram (CMD) reach- 
ing about 2 magnitudes below the oldest main sequence 
turnoff. By analysing the distribution of stars near the base 
of the red giant branch (RGB), they determined that the 
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Table 1. Observing Log 
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first generation of stars in Leo II was coeval with the for- 
mation of Galactic Globular clusters, and nearly half of the 
stars formed during a period of star formation lasting about 
4 Gyr, with the typical star forming about 9 Gyr ago. They 
found a negligible rate of star formation in the last ~ 7 Gyr. 
Subsequent re-analyses of these data confirmed this scenario 
and found a star formation history (SFH) dominated by old 
stellar populations with a low star forrnation rate in the last 
8 Gy r ^Hernandez et al.ll2000l : lDolphin|[20o3 : iDolphin et al] 
I2OO5I ) . A recent wide-field study has shown a gradient in the 
HB morphology and the mean age of stellar populations, 
with a signific ant population younger than 8 Gyr found only 
at the centre (jKomivama et aLll2007l l. 

The presence of a small intermediate age population 
in Leo II is also indicated by a smal l number of C stars 
dAaronson. Hodge, fc Olszewskil 19831: Aaronson &i Mould! 



19851: 



Azzopardi. Leaueux fc Westerlundl Il985h . 
Azzopardi et al.l (Il985| ) list six certain C stars and 



one candidate, while lAzzopardil (|2000l ) stated they found 
two new ones but without providing further details. 

As part of an imaging study of Local Group (LG) dwarf 
galaxies in the near infrared (NIR), we have undertaken a 
study of the evolved stellar populations in Leo II using wide- 
area NIR imaging. The main goals are to study the metal- 
licity distribution of red giant stars and to obtain J, H, 
and Ks magnitudes of as ymptotic giant branch (AGB) stars 
iGuUieuszik et ahllioOTal lbh. The main advantage of NIR ob- 
servations for studying AGB stars over the alternative search 
technique based on intermediate band imaging in the optical 
(I Albert. Demers fc Kunkelll2000l : iBattinelli fc Demersll2000l : 
iNowotnv et al .l200lh . is that the spectral e nergy distribution 
of coo l AGB stars peaks in the NIR (e.g.. iGullieuszik et al.l 
l2007al ). Also, bolometric corrections are smaller and more 
precise in the NIR making comparison with theoretical 
quantities easier. In addition, the (foreground and inter- 
nal) extinction is much l ower in the NIR than in the optical 
l|Rieke fc Lebofskvl[l985l ). 



2 OBSERVATIONS AND REDUCTION 

Observations of Leo II were carried out in April 19-20, 2005 
using the wide-field NIR camera WFCAM at the UKIRT 
telescope on Mauna Kea, Hawaii. WFCAM uses 4 Rock- 
well Hawaii-II HgCdTe detectors with a pixel scale of 0'.'4 
pixel" ^. The size of each array is 2048 x 2048 pixels, corre- 
sponding to about 13^6 x 13^6, and the 4 arrays are sepa- 
rated by a gap comprising 94% of the detector size. Since the 
area encompassed by one array is comparable to the galaxy 
size , we used only one array (No. 3) to observe Leo II. The 
remaining detectors provide a useful estimate of the fore- 
ground/background counts. 

Each "observation" (or pointing) consisted of 4 micro- 



stepped 10 s images (two 5 s coadds) on a 9-points jitter pat- 
tern, totalling a 360 s exposure time. In total, 6 observations 
were obtained in J and H and 10 in the Ks band, yielding 
a total on-target integration time of 36 min in J and H and 
1 h in Ka ■ The 2x2 "small" microstepping was used to ob- 
tain a better spatial sampling. The observing log is given in 
Table [U 

The raw data were processed using the WFCAM 
pipeline provided by the VISTA Data Flo w System Project , 
to which the reader is referred for details (|Dve et al.ll2006l ). 
The pipeline combines the micro-stepped images in each 
band into 4kx4k "Leavstack" oversampled images with a 
spatial resolution twice that of the original raw images, 
i.e. [/.'2 pixel"'^. The pipeline produ cts are astrometrically 
calib rated using the ZPN projection jCalabretta fc GreisenI 
2002t) and the 2MA SS Point Source Catalogue (PSC, 



Skrutskie et"al] I2OO6I ) as a reference, with a final system- 



atic accuracy of the order (/.'l. In our analysis, we made 
use only of the array where Leo II had been centred (No. 3) 
along with a second array (No. 2) used to estimate the con- 
tribution of the foreground Galactic stars and background 
galaxies. 

Point spread function (PSF) photometry was performed 
on the individual oversampled images (6 in J and H , 
10 in Ks) using the allstar/allframe (jStetsonl Il987l . 
119941 ). The PSF was generated with a Penny function with 
quadratic dependence on the position on the frame. The final 
catalogue includes instrumental PSF magnitudes of objects 
detected in at least 2 images in 2 bands. 

The positions of the sources in the raw photometric 
catalogue were converted from pixels to the J2000 equatorial 
system using the astrometric calibration provided by the 
pipeline and IRAfQ tasks including support for the ZPN 
projection. 

The final, calibrated NIR photometric catalogue of stars 
in Leo II is provided in the electronic version of the journal. 
A few lines are presented in Table [5] to illustrate its content. 

Artificial star experiments were also performed to eval- 
uate the photometric errors and completeness of our pho- 
tometry. We performed 20 test runs adding ~ 2000 stars on 
a 2fc X 2fc portion of the frames. The input magnitudes and 
colours were randomly generated to reproduce the RGB of 
Leo II. The results of our experiments (completeness factor 
and internal r.m.s. photometric errors) are shown in Fig. [21 
The completeness factor is larger than 50% for magnitudes 
brighter than Ks — 20. Note, however, that most results of 
this paper are based on photometry of stars brighter than 
Ks — 18, for which we have a completeness factor ~ 100% 
and photometric errors smaller than 0.02 mag. 

Possible spatial variations in the completeness factor 
and photometric errors were investigated by repeating the 
analysis of artificial star experiments for different regions of 
our frames, and no significant variation was found. In fact, 
our scientific photometric catalogue only contains --^ 3200 
objects, hence we expect spatially-varying crowding effects 



^ The Image Reduction and Analysis Facility (IRAF) software 
is provided by the National Optical Astronomy Observatories 
(NOAO), which is operated by the Association of Universities 
for Research in Astronomy (AURA), Inc., under contract to the 
National Science Foundation. 
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Figure 1. Comparison of Leo II WFCAM instrumental magni- 
tudes corrected using the [d^ et al. colour terms with 
2MASS photometry, for stars in common with the 2MASS PSC. 
No residual colour terms are detectable. In each panel the median 
difference and the standard deviation of the data are shown. 
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Figure 2. Completeness (lower panel) and photometric errors 
(upper panel) in our Leo II images as derived from artificial star 
experiments, plotted against the input magnitudes of simulated 
stars. 



Table 2. The NIR catalogue of Leo II stars over WFCAM ar- 
ray No. 3. A few lines are shown here for guidance regarding its 
form and content, while the full catalogue is available from the 
electronic edition of the journal. . 
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to be unimportant. The effect of a photometric bias to- 
wards bright er retrieved magnitu des, caused by photometric 
blends fe.g.. iGallart et al.lll996l l. was also investigated and 
found negligible (in our experiments the mean difference is 
always less than 0.01 mag for stars brighter than Ks ~ 18). 



3 COLOUR-MAGNITUDE DIAGRAMS 

Figure [3] presents our J — Ks, Ks colour- magnitude diagram 
of Leo II, along with optical-NIR CMDs and an optical 
B — V,V CMD. A selection based on the SHARP parameter 
was applied to remove noise peaks, diffuse objects, and other 
spurious detections . The optical-NIR CMDs and the opti- 
cal CMD were obtained by adding data ob tained with the 
EMM I camera at the NTT at ESO/La Silla (|Momanv et al.1 
I2OO7I . Rizzi et al., in prep.). Since the optical photometry 
refers to a smaller central area (9^1 x 9!l in V , 6.'2 x 6^2 in 
B), all diagrams involving optical data refer to a subset of 
our WFCAM catalogue. 



Our NIR CMD of Leo II shows a well populated RGB, 
with an T RGB clearly visible at Ks ~ 16. Two 10 Gyr old 
isochrones (|Girardi et al. | |2002| ). with metallicity Z = 0.001 
and Z = 0.0004 (corresponding to [Fe/H] ~ -1.3 and [Fe/H] 
~ —1.7), are superimposed to the CMD. The CMD almost 
reaches the level of the HE, identified with the tail of faint 
stars with colours bluer than the RGB, at a magnitude 
which is consistent w ith Vhb = 22.18 ±0.18 measured by 
iMigheU fc Richl (f 19961 ). However, these stars will not be fur- 
ther analysed since they fall close to the detection limit in 
the CMD. 

The sequence of stars brighter than the TRGB (Ks ~ 
15) is identified with upper-AGB stars belonging to an 
intermediate-age stellar component. These stars appear 
to coincide with the sample of C stars identified by 
lAzzopardi et al.l (| 19851 ). Note that they are more luminous 
than the TRGB only in the J — Ks, Ks diagram, while they 
appear progressively fainter in bluer photometric bands. 
This C star population is quite small, in agreement with 
the low star formation rate of LeoH at intermediate ages. 
The AGE population will be further discussed in Sects. [6] 
and [71 The cand i date C star with uncertain classification in 
lAzzopardi et al.l (| 19851 ) does not share the spectral energy 
distribution of the C stars in the different CMDs. Visual 
inspection of our images confirms that it as a background 
galaxy. 

Note that both the prominent RGB sequence and the 
AGE component are absent in the CMD of an outer field of 
Leo II, obtained from the detector No. 2 of WFCAM (Fig.g}. 
This field is located at about 26' from the centre of Leo II, 
corresponding to about 3 tidal radii. A comparison with a 
simulation of the Galactic foreground, obtained for the sam e 
field-of-view using the trilegal code (|Girardi et al.ll2005l ). 
is shown in the right panel of Fig. [4] The comparison in- 
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Figure 3. Colour- magnitude diagrams of Leo II from NIR and combined optical-NIR photometry. In all diagrams, C stars spectroscop- 
ically identified by lAzzopardi et all l|l985l ) are shown as filled circles, while their C s t ar un certain candidate is marked by a triangle. 
Superimposed on the J — Ks , Kg CMD are theoretical isochrones from iGirardi et al. I ||2002|'I . for an age 10 Gyr and two metallicities 
close to the metallicity of Leo II, Z = 0.001 and Z = 0.0004. In all ex cept the J — Ks, Ks CMD, we also show as asterisks the predicted 
Galactic foreground stars towards Leo II from a TRILEGAL simulation llGirardi et al.|[2005h . The He-burning stars on the "yellow plume", 
selected from the optical B — V,V diagram, are shown as open squares. They are best seen in the B — Ks, B diagram, whereas the same 
stars arc confused with the RGB in the J — Ks, Ks CMD. 



dicates that the CMD of the outer region mostly contains 
Milky Way stars (the vertical sequence with J — Ks ~ 0.8 
and bluer objects) and background galaxies (concentrated 
at J-Ks >1, Ks> 18). 

The wide baseline of optical-NIR diagrams provides the 
best separation of some CMD features. In Fig. |3] we plot 
as open squares the stars on the "yellow plume" or "verti- 
cal clump" (VC), just above the red clump, according to a 
magnitude and colour selection from the optical B — V, 1/ di- 
agram. These stars are best seen in the B — Ks , B diagram as 
a vertical sequence originating from B~Ks ~ 2.5, B ~ 22.5 



and extending up to B ~ 21; while they are hardly detected 
in the J — Ks, Ks NIR diagram, an ambiguity that is ex- 
plained by a combination of increasing photometric errors 
at faint magnitudes and an intrinsically narrow baseline. 

The distribution of these "vertical clump" stars can be 
compared in Fig.[3]with the distribution of Milky Way fore- 
ground stars in the direction of L eo II, obtained from the 
TRILEGAL code (jCirardi et al.ll2005h . The projected Galactic 
contamination (plotted as starred symbols) is insignificant 
in the yellow plume region of the B — Ks, B diagram, indi- 
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Table 3. Observed magnitude of the TRGB and distance moduli 
derived for Leo II from JHKs photometry. 



Figure 4. Left: colour-magnitude diagram of an outer field of 
Leo II, obtained from the detector No. 2 of WFCAM, pointing at 
about 26' from the centre of Leo II. The vertical sequence having 
J — Ks < 0.8 is consistent with the expected Galactic foreground, 
while the group of objects redder than J — Kg = 1 and fainter than 
Ks = 18 are unresolved background galaxies. Right panel: Milky 
Way stars i n the Leo II line of s ight as derived from a TRILEGAL 
simulation l lGirardi et al.ll2005h . for the observed field-of-vicw, 
without photometric errors. 



eating that the vertieal sequence is certainly a Leo II stellar 
population. 

In dwarf galaxies, this feature is usually attributed to 
a population of core He burning stars a few hundred Myr 
to ~ 1 Gyr old, which are the descendants of stars lo- 
cated above the old main-sequence turnoff , the so-called blue 
plume. Thus, the detection of VC stars is generally inter- 
preted as evidence of recent star formation in these galax- 
ies (e .g., in Draco: lAparicio. Carrera fc Martmez-Delgadol 
I2OOII ). This is probably the case for Leo II, where there 
is evidence for a n increasing number o f young stars to- 
ward the centre (|Komivama et ahl l2007f ). However, there 
is some evidence that the detection of a VC sequence in 
dwarf spheroidal galaxies may not be suffi cient to estab- 
lish the presence of rece nt star formation (|Momanv et al.l 
I2OO7I : iMapelh et al.|[2007[ ). since stars brighter than the HB 
ha ve been detected in globular clusters (see the case of M80 
m iFerraro et "aLlll999l ). The agreement of the "blue plume" 
frequency in 7 dwarf spheroidal galaxies (including Leo II) 
might be a hint that this population is partly comprised of 
"blue stragglers" , of which the VC population may represent 
the evolved counterparts. 



4 DISTANCE FROM THE RGB TIP 

The luminosity of the TRGB in the I band (where the de- 
pendence on the age and chemical composition of the stel- 
lar population is at a mini mum) has long been used as 
a valuable standard candle (|Da Costa fc Armandrof3ll990l : 
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iLee et all Il993h . In the NIR, the TRGB luminosity de- 
pends on age and metallicity in a more complex way. 
For instance, the Ks-h&nd magnitude of intermediate-age 
stars at the TRGB is fainter than that of old stars; while 
the TRG B Ks luminosity rises with increasing metallic- 
ity (e.g., ISalaris fc Girardil |2005| ). For a population that 
becomes more metal-rich with time as a result of galaxy 
chemical evolution, the two effects can partly balance. In 
iGuUieuszik et al.1 (|2007al 'l we showed that, if the galaxy's 
SFH can be (even roughly) estimated, quite accurate dis- 
tance determinations can be obtained from the TRGB at 
NIR wavelengths. 

The method is applied here to measure the distance to 
Leo II independently of optical measurements. This will also 
be a useful test of the reliability of distance estimates in the 
NIR domain, which is important for next-generation instru- 
ments operating mainly at NIR wavelengths (e.g., JWST, 
adaptive optics at Extremely Large Telescopes). 

We estimated the magnitude of the TRGB of Leo II 
by fitting its Ks-hand luminosity function to a step func- 
tion convolved with a Gaussian kernel representative of the 
photometri c errors. This method , extensively applied by our 
group fe.g.. lMomanv et al.|[20o3 ). was found to give consis- 
tent resu lts within la w i th the Max imum Likelihood Algo- 
rithm of iMakarov et al.l (|2006h (see lRizzi eraLll2007l ). The 
resulting J, H, and Ks magnitudes of the TRGB are given in 
TableO The errors associated to these magnitudes are dom- 
inated by the uncertainty on the absolute photometric cali- 
bration, because the error resulting from the TRGB fitting 
algorithm is less than 0.01 mag, and the internal photometric 
errors at the level of the TRGB are negligible (see Sect. [2]). 
The observed magnitude were cor rected for extinction using 
a reddening Eb-v = 0.03 and the lRieke fc Lebofskvl (|l985i ) 
reddening law. 

To derive the distance to Leo II, the JHKs TRGB mag- 
nitudes were compared with the empirical calibrations of 
j^^TRGB ^ a function of [M/H] based on Galactic globu- 
lar clusters (jValenti et al.ll2"004bl ). whose intrinsic system- 
atic error is cr = 0.16 mag. The adopte d mean metallicity 
was [M/H] = -1.73, in agreement with iKoch et al.1 (|2007h 
spectroscopy. Given the relatively old age distribution of 
Le o II, we found the popul ation corrections (calculated as 
in iGuUieuszik et al.ll2007ah to the TRGB magnitude to be 
negligible. 

The distances derived from the J, H, and Ks bands are 
also given in Table |3] The weighted mean is (m — M)o = 
21.68 ± 0.11, and intermediate be tween the "shor t " dis - 
tance (m - M)o = 21.55 obtained bv lMigheU fc RichI (119961 ) 
from the V magnitude of the HB a nd the "long" d i stance 
(m-Af)o = 21.84 ± 0.13 derived bv lBeUazzini et al.1 (|2005l ) 
from the I magnitude of the TRGB. Therefore, our determi- 



6 GuUieuszik et al. 




Figure 5. Upper panel: the photometric MetaUicity Distribu- 
tion Function of RGB stars in Leo II, constructed using global 
metalHcities [M/H] derived from {V — Kg) colours. The solid line 
is a Gaussian fit to the data, with mean [M/H] = —1.74 and 
"'[M/H] = 0.20 dex. The measurement scatter is represented by 
a Gaussian with (Tinstr = 0.06 dex {da shed curve). Lower panel: 
the spectr oscopically d e rived MDFs of iKoch et al ] l|2007h (solid 
line) and iBosle r et al.l ||2007|) (dashed histogram), both on the 
ICarretta fc GrattonI lligOTl t "scale. 



nation based on NIR magnitudes of the TRGB appears to 
be consistent, within the errors, with the optical estimates. 
The techniques explored in this paper will be useful to mea- 
sure the distance of stellar systems for which photometry of 
resolved stars will become available only in the NIR. 



5 METALLICITY 

5.1 MetaUicity distribution of RGB stars 

Given the uncertainties on the metallicity of Leo II and the 
importance of the metallicity distribution function (MDF) 
of RGB stars as a constraint on models of chemical evo- 
lution of dwarf galaxies, we used all the information from 
optical-NIR colours to investigate the metallicities of stars 
in Leo II. This parameter is the most appropriate to esti- 
mate the metallicities of dwarf spheroidal galaxies (having 
[a/Fe] ratios close to solar) by comparison with the photo- 
metric properties of Milky Way globular clusters, which gen- 
erally show an overabundance of a-elements r elative to iron 
that is a function of the cluster metallicity (see lGeisler et all 
I2OO7I . and refs. therein). 

The photometric MDF obtained for red giant stars in 
Leo II down to 2 mag below the TRGB, is shown in Fig. [5] 
(upper panel). The distribution is well described by a Gaus- 
sian function centred at [M/H] = —1.74 and with a measured 
dispersion of 0.20 dex. The internal error in the same mag- 



nitude range was evaluated by applying the same method 
to a synthetic CMD simulating a thin RGB, taking into ac- 
count the results of artificial star experiments. The recovered 
metallicities have a Gaussian distribution with a dispersion 
0.06 mag, assumed to be representative of the internal error 
of our metallicity measurements. By quadratically subtract- 
ing this internal error from the measured width of the MDF, 
we obtain a corrected dispersion 0.19 dex for the photomet- 
ric MDF of Leo II. 

Our photometric MDF shown in Fig.[5]is representative 
of the true MDF only for stars as old as the Galactic GCs 
(~ 1 2.5 Gyr). However, a typical Leo II star is 9 Gyr old 
(e.g. iMighell fc Rich|[l996l '), hence slightly bluer than glob- 
ular cluster stars of the same metallicity. Therefore, as in 
I GuUieuszik et all (j2007al l. we used theoretical isochrones to 
construct contours of constant V — Ks colours of RGB stars, 
as a function of both stellar age and metallicity, and correct 
the measured metallicity for the age effect. This is done by 
estimating the metallicity of a 9 Gyr old star having the 
same colour as a 12.5 Gyr star with [M/H] = —1.74, i.e. 
the mean of the MDF. This differential approach overcomes 
any possible problems with the absolute calibration of the 
isochrone colours. With this assumption, the age correction 
results to be A[M/H] = 0.10, which is very small and com- 
parable with the absolute uncertainty of our method. Our 
major simplification is that all stars have the same age, but 
it is adequate to calculate the mean metallicity of Leo II 
stars. By applying this correction, the mean metallicity of 
Leo II turns out to be [M/H] = -1.64±0.06 (random) ±0.17 
(systematic) . 

5.2 Comparison with spectroscopy 

The derived metallicity is in excellent agreement with the 
two rece nt spectroscopic resul ts [Fe/H] = —1.73 and = 
-1.59 bv lKoch et all (|2007D and lBosler etai] (|2007l ). respec- 
tively. Our metallicity distribution of Leo II RGB stars, as 
inferred from V — Ks colours, is co mpared with the spe ctro- 
scopic meta l licity distributions from lBosler et ahl (|2007h and 
iKoch et all |20^ in Fig. [S] (lower panel). The distributions 
are basically consistent, except for a slightly lower mean 
metal abundance from photometry. In particular, the range 
in metallicity (FWHM of the distributions) is comparable. 

The significant overlap bet ween our sample of RGB 
stars and the catalogues of iBosler et all |20o3) and 
iKoch et all (|2007l ') allows us a direct comparison of photo- 
metric and spectroscopic metallicity estimates on a star-by- 
star basis. 

Figure [6| sho ws a comparison with the results of 
IBosler et aL^(f2007^ ■ for 71 stars in common with our sam- 
ple. In the metallicity range typical of Leo II stars, the 
two calibrations adopted by those authors (as a function 
of [Fe/H] and [Ca/H]) yield [Fe/H] ~ [Ca/H]. Indeed, the 
relations presented in Fig.|S]for the two scales are quite sim- 
ilar. For both scales, the overall agreement between spec- 
troscopic and photometric metallicities appears to be good. 
The most noteworthy difference is for the bluest stars, whose 
spread is higher and spectroscopic metallicities are system- 
atically higher. A possible explanation is that the photo- 
metric metallicities are underestimated because some stars 
are actually younger (hence bluer) than old RGB stars. One 
alternative possibility is that calibration uncertainties and 
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Figure 6. Comparison of our me asured photometric metallicities 
with the spectroscopic results of iBosler et al.l l l2007f) . In the up- 
per panel we plot the spectroscopic data on the [Ca/H] scale of 
the authors, while the calibration based on [Fe/H] is shown in the 
lower panel. The error crosses represent the mean uncertainties of 
spectroscopic metallicities, and the Icr internal errors of our pho- 
tometric determinations (not including age effects and systematic 
uncertainties). 



internal errors aflect in some way the spectroscopic measure- 
ments at low metallicity. 

A comparison of our photometric metallicities (with 
no age corr e ction applied) with the spectroscopic results of 
iKoch et all (|2007l ') is presented in Fig. [T] The 41 stars in 
common with our sample are divid ed in 3 age interva ls, us- 
ing the age estimates publis hed by Koch et al.l (|2007h . The 
stars with ages (as derived bv lKoch et al.ll2007 ) greater than 
7 Gyr, and most of the stars with ages in the range 4 to 7 
Gyr, are broadly consistent, with a large scatter, with the 
bisector in Fig. [T] i.e. compatible with the age of Galactic 
globular clusters. 

We have also plotted the age-corrected relations (repre- 
sented by different lines in Fig. [7| , by calculating the metal- 
licity shifts to be applied to our photometric measures for 
young stellar populations. Assuming an age of 9, 7, and 4 
Gyr, the corrections are A[M/H] = 0.10, 0.20, and 0.41, re- 
spectively. Indeed, the location of stars younger than 4 Gyr 
seems consistent with the expected relation for 4 Gyr old 
stars. In general, however, we note a sizeable scatter, even 
co nsidering stars with in each age bin, and the ages estimated 
bv lKoch et all (|2007h do not appear to be closely correlated 
with the age-corrected relations in Fig. [7 Overall, th e num - 



ber of young stars in Leo II derived by iKoch et al.l (|2007l ) 
appears to be larger than su ggested by SFH recons tructions 
based on HST photometry (|Hernandez et al..,2000i : iDolphirJ 
|2002| . Rizzi et al., in prep.), and the mean age of RGB stars 
is younger. We note that our age corrections, based on the 
larger baseline of V — Ks colours, may give more p r ecise age 
ranking than the g — i colour used by iKoch et al.l (|2007l ) . 



Figure 7. Com parison of individ ual stellar metallicities with 
those derived by iKoch et al. using different symbols for 

stars in different age intervals, as estimated by Koch and coll. 
The lines represent the expected relations when the effects of age 
on V — Kg colours are taken into account. For example, 4 Gyr old 
stars arc expected to follow the long-dashed line. Error bars as in 
Fig.E] 



A direct comparison of spectroscopic metallicities, 
which would be interesting to assess the accuracy of spec- 
troscopic metallicities, is not possi ble due to t he ab sence 
of overlap between the samples of iKoch et al.l (|2007l ) and 
iBosler et all ((2003). 



6 TWO-COLOUR DIAGRAMS: SELECTION 
OF AGB STARS 

In this section we present the two-colours diagram, which 
is used to select AGB stars in Leo II. This diagram is a 
powerful tool to separate the foreground Milky Way stellar 
populatio n, and allows a separation of carbon and oxygen - 
rich stars (|Aaronson fc Mould|[l985l : iBesseU fc Brettlll988l ). 

Figure [5] shows the NIR two-colour diagrams of Leo II 
and the external field. We selected only stars brighter than 
the TRGB {K^ = 15.75) to exclude RGB stars. Stars are lo- 
cated in well defined sequences, within the regions outlined 
in Fig. [8] All stars located in re gions 2 and 3 are found along 
the dwarf stars locus defined by Bessell fc BrettI d 19881 ) . The 



number of stars in each region are also given in Fig.|8]for the 
field centred on Leo II and the external field. The number of 
stars in regions 2 and 3 are the same, within statistical fiuc- 
tuations, in the two fields. We can therefore conclude that all 
stars in regions 2 and 3 are Milky Way dwarfs. On the other 
hand, stars in region 1 and 4 are found only in the field cen- 
tred on Leo II, and we can conclude that they are all Leo II 
members. Being brighter than the TRGB, they can only 
be AGB stars. For stellar populations younger than those 
present in Leo II, core He-burning red supergiants should be 
also considered. 
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Figure 8. Left Panel: the two-colour diagram of Leo II stars brighter than the TRGB {Ks = 15.75), with superimposed the regions we 
have used to discriminate stars in Leo II from those in the Milky Way: region (1) are probable Leo II O-rich AGB stars, regions (2) and 
(3) are dwarf Galactic stars, region (4) is populat ed by C stars in Leo I I. Different symbols indicate stars in different regions. The loci of 
giant stars and main-sequence dwarf stars (from iBessell fc Brett|[l988l ) are shown as a dashed and solid line, respectively. Right Panel: 
the same, for the outer field. Note the absence of stars belonging to Leo II. 



All of the five stars in reg ion 4 were identified as C 
stars by lAzzopardi et al.l (|l985 ). One more C star in their 
catalogue (ALW5) is slightly fainter than the TRGB (see 
Fig. 13]) and hence is not marked in Fig. [51 If plotted, it 
would fall in region 3. 

Figure [5] shows the location in the CMD of the stars 
classified using the NIR two-colour diagram. Stars in the re- 
gions 1 and 4 are consistent with the expected loci of M and 
C stars respectively, as judged from the location of AGB 
stars with spe ctroscopic classification i n the CMD of For- 
nax dSph f see iGuUieuszik et al.l [2[)07al . and ref.s therein). 
The C star population, in particular, agrees well with the 
mean colour-magni tude relation for C stars in LG dwarf 
galaxies derived by iTotten et al.l (|2000l ). scaled to the dis- 
tance of Leo II discu ssed in Sect. 31 N I R ph otometry of all 
C stars classified by lAzzopardi et al.l (|l985l ) and probable 
O-rich AGB stars selected by us in region 1, is given in Ta- 
bled 

In the following , we c onsider all the C stars identified 
by lAzzopardi et"aL (|l985l ). including star ALW5, which is 
fainter than the TRGB and was therefore not included in 
Fig. |8] The objects ALW2, misidentified by lAzzopardi et"al] 
(|l985h as a C star, is not included in our analysis. All the 
remaining stars by lAzzopardi et al.l (|l985l ) are compatible 
with our C star selection and our observations cover all 
Leo II, nearly out to the tidal radius. We therefore conclude 
that the complete population of Leo II C stars in the central 
13^6 X 13.'6 area covered by our observations, is f ormed by 6 
object s (excluding ALW2). We finally note that lAzzopardil 
(|2000l ) stated they found 2 new ones but without providing 
further details. In our selection we have no indications for 
the presence of other objects in addition to the 7 discussed 



here. Finally we consider as O-rich AGB stars the 7 stars 
found in region 1 in Fig. [8] 



7 COMPARISON WITH THEORETICAL 
MODELS 

A distinctive feature of present observations is that they 
sample quite completely the optically-visible AGB popula- 
tion of Leo II in the surveyed area. The only AGB stars 
expected not to be present in our data are those so strongly 
absorbed by circumstellar dust to become invisible even in 
the NIR. 

Such complete catalogues of AGB stars in nearby 
galaxies are rare. The best such data are no doubt 
those for the LMC and SMC, fully sampled in the 
IJHK^ bands of DENIS an d 2MASS (see lCioni et"aLlll999l : 
iNikolaev fc Weinberg! bOOOl'l an d now being sample d in the 
mid-IR (e.g iBlum et aL niooel : iBolatto et al.l [2OO7I) . Com- 
pared to these galaxies, Leo II is more metal-poor (Sect. [5]), 
and presents a much simpler history of star formation, con- 
centrated at old ages. These particularities provide us with 
a unique opportunity to test present-day AGB models in the 
interval of low masses and low metallicities. 



7.1 Simulating the photometry 

Having this goal in mind, we will try to fit the Leo II 
observed AGB population with the recent set of ther- 
mally pulsing AGB (TP -AGB) evolutionar y tracks from 
iMarigo fc Girardil (|2007h . Added to the ICirardi et al] 
( 2OOOI ') tracks for the pre- TP- AGB evolution, they 
are converted to stellar isochrones as described in 
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Figure 9. The position of Leo II and foreground stars in the CMD according to the classification shown in Fig. [8] for stars brighter than 
the TRGB. Left: stars in the regions 1 {squares) and 4 (filled circles). The open t riangle represents ALW 2, a misidentified background 
galaxy, while the filled triangle is a C star found just below the TR GB (ALW5, Azzopardi et al1ll985^ . The solid curve is the mean 
colour-magnitude relation for C star in Local Group dwarf galaxies llTotten et al. scaled to the distance of Leo II. Right panel: 

the same, for probable Milky Way stars in the regions 2 (open squares) and 3 (open circles) of the two-colour diagram. 



Table 4. NIR photometry of Leo II A GB stars. The identifie rs are those in our photometric 
catalogue. For the C stars, the names in lAzzopardi et al.l (Il985l ) are also given. 



ID 


a (J2000) 


S (J2000) 


J 


H 


Ks 


type 


note 


1661 


11:13:12.82 


4-22:11:14.1 


16.231 


15.284 


14.788 


C 


ALWl 


1671 


11:13:20.64 


4-22:11:16.3 


16.366 


15.351 


14.812 


C 


ALW3 


659 


11:13:23.48 


-1-22:07:58.4 


16.542 


15.674 


15.308 


C 


ALW4 


828 


11:13:23.97 


H-22:08:29.3 


17.054 


16.365 


16.122 


C 


ALW5 


1089 


11:13:29.39 


4-22:09:14.2 


15.963 


15.114 


14.876 


C 


ALW6 


1032 


11:13:31.78 


4-22:09:06.1 


16.456 


15.679 


15.433 


C 


ALW7 


781 


11:13:20.83 


4-22:08:22.9 


16.005 


15.184 


15.027 







1215 


11:13:35.81 


4-22:09:35.1 


16.252 


15.471 


15.351 







1509 


11:13:29.24 


4-22:10:32.9 


16.346 


15.531 


15.375 







1597 


11:13:23.16 


4-22:10:56.0 


16.487 


15.740 


15.574 







1873 


11:13:53.43 


-^22:12:43.5 


16.229 


15.510 


15.404 







1877 


11:13:52.77 


4-22:12:45.4 


15.855 


15.113 


14.989 







1904 


11:13:29.17 


4-22:13:02.7 


16.189 


15.387 


15.236 








iMarigo et all (|2008l ') and fed to the trilegal pop- 
ulation synthesis code for simula ting the photome- 
try of resolved stellar populations (|Girardi et al.l |2005| . 
|http : // trilegal . kuleuven . be/ ) . Since the details of the 
TP- AGB imple mentation in trilegal a re provided in sepa- 
rate papers (e.g. lGirardi fc Marigd[2007l , and work in prepa- 
ration), suffice here to recall the basic aspects of the model 
simulations: 



• The Mi l ky W ay foreground is simulated as in 
iGirardi et al.l |200^, including the main disk and halo com- 
ponents and for the same area of our observations. 



• The Leo II galaxy is set at a distance of 205 Kpc (this 
paper). Reddening is ignored since it is negligible. 

• We fix the size of the simulations by reproducing the 
star counts in the upper 2 magnitudes of the RGB - for 
which our observations are quite complete - together with 
the relative star formation rate (SFR). 

• The relative SFR of Leo II is taken from two differ- 
ent sources, a s dep icted in the upper panel of Fig. 1101 
iDolphin et"al] (120051 ) and Rizzi et al. (in prep.). In both 
cases the SFR is derived from the inversion of a deep CMD 
from HST. Although the error bars in these SFR determi- 
nations are quite significant, these SFRs show that the bulk 
of star formation in Leo II was confined to ages larger than 
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2 Gyr, and concentrated at > 4 Gyr. More details are given 
later. 

• The metallicity is taken from the photometric determi- 
nation of Sect. [21 i.e. we use a Gaussian distribution of mean 
[M/H]= —1.64 and dispersion 0.2 dex. This metallicity dis- 
tribution is assumed to be the same for all ages, which is 
likely a good a pproximation since the age-met allicity rela- 
tion derived bv lKoch et al.l (|2007l . their figure 11) is essen- 
tially flat for ages above 5 Gyr. 

• Each simulated star is conver ted in the 2MASS system 
using an updated version of the iBonatto. Bica fc Girardil 
l|2004h transformations. In particular, t he transformations 
for carbon stars are now derived from Loidl et all (|200ll ) 
spectra. In addition, using iGroenewegenl (|2006l ) tables we 
correct the photometry for the effect of circumstellar dust 
in mass-losing AGE stars. The 60% Silicate -I- 40% AlOx 
and 85% AMC -I- 15% SiC dust mixtures are assumed for 
O-rich and C-rich stars, respectively. 

• For TP-AGB stars, the pulse cycle luminosity and T^s 
variations are also simulated; long period variability is not. 
For the high-amplitude Miras, variability may provide an 
addit ional scatter of about 1 mag (see, e.g., ICioni et al.l 
I2OO3I ) in the K band. Our simulations however predict that 
most of the AGB stars are first overtone pulsators, which 
have much smaller amplitudes. 

• Photometric errors and completeness are simulated us- 
ing the relations derived in Fig. [2] 

• To reduce the statistic fluctuations in the numbers of 
predicted stars, each simulation is run at least 50 times with 
different random seeds. When we refer to the "expected 
numbers" of each kind of star, we are actually referring 
to the mean values and standard deviations obtained from 
these many runs. 
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Figu re 10. Top panel: Th e star-formation histories used in this 
work: lDoIphin et al.l bOOSl . dotted line) and Rizzi et al. (in prep.) 
{solid line) as a function of age. Bottom panel: The production 
(number of stars per time interval) of different types of stars in 
our models as a function of age, for a galaxy model forming stars 
at a constant rate (in mass per unit time) from to 15 Gyr, at 
a constant Z = 0.0004 metallicity. The stellar kinds plotted are 
RGB stars within 2 mag of the TRGB (multiplied by 0.4; dotted 
line), and both O-rich giant stars above the TRGB {dashed line) 
and C-rich giants {solid line). 



Our simulations also take into account the error bars in 
iDolphin et al.l (|2005h and Rizzi et al. (in prep.) determina- 
tions of the SFH. These error bars reflect both the intrinsic 
errors in the method of SFR-recovery, and the small number 
statistics of the original HST data from which the SFR is 
derived. For each age interval, we use a random SFR value 
drawn from a normal distribution centred at the mean SFR 
and with the appropriate value of a. Since for some age bins 
the 1(7 error bars are comparable to the mean SFR, the neg- 
ative values of SFR obtained at the youngest age bins are set 
to zero. This produces a distribution of SFR values that is 
non-symmetrical around the mean values, especially at the 
youngest age bins. 

One of such simulations is shown in Fig. 1111 which re- 
sembles very much Fig. [3l First, a key to unde rstand the sim- 
ulations is given by the lMarigo et all (|2008l ) isochrones for 
a few selected ages. They show the location of the TP-AGB 
phase at quiescent phases of H-shell burning, where these 
stars spend about 70 % of their life. The remaining ~ 30 % 
is spent at phases of lower luminosity (and higher Tes) after 
the occurrence of He-shell pulses. The result is that the TP- 
AGB stars in the simulation are typically found above the 
TRGB, in the same region defined by the isochrones, but a 
tail of such objects (both C- and O-rich) extends down to 
almost 2 mag below it. Moreover, a significant fraction of 
the O-rich giants found above the TRGB are not genuine 
TP-AGB stars, but either early-AGB stars or (more rarely. 



in the case of Leo II) core He burning stars belonging to the 
youngest populations. 

Comparing the simulated data points with the observed 
ones of Fig. [31 one notices that the main stellar features are 
accounted for by the model. The simulations do not contain 
the objects observed at the bottom right part of the dia- 
gram {Ks > 18, J — Ks > 1.2), which likely corresp ond to 
background galaxies (e.g. iNikolaev fc Weinber3l2000l ). Sim- 
ulated carbon stars are on average bluer than the observed 
ones, although a few very red dust-enshrouded objects are 
present in our simulations. There are also minor offsets in 
colours of O-rich stars, that cannot be appreciated in Fig. [TT] 
because they are of the order of 0.05 mag. It is interesting 
to note that 3 out of the 5 stars expected to be C-rich, are 
located along the O-rich sequence, i.e., not exhibiting red- 
der J — Ks colours than those of the O-rich stars with the 
same luminosity. This prediction is consistent with the ob- 
served colour-magnitude diagram presented in Fig. [3l (top- 
left panel) , where 3 C-rich objects are just seen to lie on the 
O-rich giant branch. On the theoretical ground this feature 
is explained considering that at low metallicity the cooling 
effect on the stellar atmosphere due to the carbon-enhanced 
molecular opacity becomes less efficient, as illustrated by 
iMarigo fc Girardil (|2007l . their figure 7). 
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Figure 11. An exam ple of simulated CMDs for Leo II, using the 
iDolphin et al. SFH. In the electronic version of this paper, 

different colours mark different kinds of stars, namely: Milky Way 
disk (green crosses) and halo (magenta dots), Leo II pre-TP-AGB 
stars (dark dots), Leo II early-AGB stars above the TRGB mag- 
nitude (cyan crosses), and Leo II TP- AGE stars bot h O-rich (blue 
square s) and C-rich (red circles) . We also plot the iMarigo et alj 
1I2OO8I ') isochrones for Z = 0.0004 and ages of 1, 2, 4, and 8 Gyr, 
shifted by the Leo II distance modulus. Dark lines mark phases 
previous to the TP-AGB, blue lines the O-rich TP- AGS, and red 
lines the C-rich one. The TP-AGB lines correspond to quiescent 
phases of evolution. 

7.2 Comparing foreground and RGB counts 

We find that the expected number of foreground stars in 
our 0.052 deg'^ area, limited to the 13 < Ka < 16 magnitude 
interval, is 38.6 ± 6.8; this is well compatible (within la) 
with the 44 objects observed in regions 2 and 3 of Fig. [8] at 
Ks < 16. This agreement is just expected, since one of the 
deep fields used to calibra te trilegal - namely the CDFS 
l|Groenewegen et al.ll20o3 ) - is located at the same galactic 
longitude and at a similar latitude from the Galactic Plane 
(i.e. i = 220?0, b = -53?9) as Leolfl. Therefore, we would 
expect that the typical errors in the predicted number counts 
at the position of Leo II {£ = 220?2, b = -f67?2) are similar 
to those of t he CDFS, i.e. of just ~ 10 % down to Ks ~ 18 
(see Fig. 6 in lGirardi et~aLll2[)05l ). 

As the simulations predict the correct number of fore- 
ground stars at 13 < Ks < 16, they can also be used to 
infer the field contamination of other CMD regions. We find 
that a total of 62.0 ± 7.2 foreground stars are expected to 
contaminate the uppermost 2 mag of the RGB. The total ob- 
served number is 254. Therefore, Leo II genuine RGB stars 
are expected to be about 192 and outnumber the foreground 
contaminants in the upper part of the RGB by a factor of 
about 4. If we assume the possible errors in the number of 

^ Note that number counts in NIR bands are very much symmet- 
rical with respect to the Galactic Plane, at least for \b\ > 20 deg. 



simulated foreground stars to be of the order of 20 percent, 
they will have just modest consequences (errors of the order 
of 8 %) in determining the total Leo II mass to be simulated. 



7.3 Comparing AGB counts 

We proceed with Leo II simulations with a total mass scaled 
such as that 192 ± 14 RGB stars are produced in the upper 
2 magnitudes of the RGB. The results concerning the AGB 
are as follows. 

For the Rizzi et al. (in prep.) SFR, we expect to find 
43.9 ± 6.6 O-rich giants above the TRGB (with 29.4 ± 5.5 
being genuine TP-AGB stars). The C-rich AGB stars are 
4.4 ±2.4 (with 3. 4 ±2.2 above the TP-AGB). In comparison, 
the Leo II data presents 7 and 6 of such stars, respectively. 
There is a clear excess of O-rich giants in the simulations, 
by a factor of about 6, which is extremely unlikely to be due 
to statistical fluctuations. For the C stars, instead, getting 
6 stars out of an expected number of 4.4 is well inside the 
67 % confiden ce level (CL) o f a Poi sson distribution. 

Using the lDolphin et all (|2005l ) SFH, 42.0 ± 5.5 O-rich 
giants are predicted above the TRGB (26.4 ± 4.7 genuine 
O-rich TP-AGB ones), and 8.3 ±2.7 C-rich (5.7 ±2.1 above 
the TAGB). The excess of O-rich giants is again of a factor 
of about 6. For the C stars, the 8.3 predicted stars are again 
inside the 67 % CL of a Poisson distribution of the 6 observed 
ones. Therefore, also in this case the observed C stars are 
compatible with the model predictions. 

To understand why models using Rizzi's et al. SFH 
present about half of the C stars as compared to the 
iDolphin et al] (|2005h case, it is instructive to compare the 
two panels of Fig. 1101 the top panel showing the relative 
SFRs and the bottom one showing the age distribution of 
different stars for a model galaxy of the same metallicity but 
with constant SFR. The bottom panel shows that the max- 
imum age for the formation of carbon stars, tc™i is close to 
6 Gyr ago. This limit is actually determi ned by the lifetime 
of the least massive TP-AGB stars in the lMarigo fc Girardil 
(|2007. ') models to experience the third dredge up events, with 
1.0 Mq (6 Gyr). Below tQ^^ the C stars predominate, above 
it they are simply absent and the relative number of O-rich 
TP-AGB stars increases. A substantial fraction of the star 
formation in Leo II has occurred close to t^^^, and this deter- 
mines a marked dependence of the C star co unts on the de- 
tails o f the SFR at this age interval. Since the lDolphin et al] 
(|2005l ) SFR presents a marked episode of star formation be- 
tween 6 and 8 Gyr, this determines the large number of C 
stars of the corresponding model. 

Were the mass limit for the dredge-up to occur just 
10 % (or 0.1 Mq) different, t^"^^ would change by as much 
as ~ 40 %. This would impact very much on the predicted 
numbers of C-type stars. Also the numbers of O-rich AGB 
stars would be affected, although somewhat less, since they 
reflect the complete SFR up to ages of 15 Gyr. The simple 
fact that the numbers of predicted C-type AGB stars turn 
out to be consistent with observations within the 95 % CL 
of a Poisson distribution, would be indicating that the min- 
imum mass for the formation of C stars at low metallicities 
is indeed close to 1.0 Mq. This is an important indication 
for the theoretical modeling of AGB stars. We recall that at 
LMC metallicities, the same mass limit is closer to 1.4 Mq. 
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Figure 12. Simulated LFs for luminous stars in Leo II, separated 
as O-rich giants above Kb = 15.75 {top panel) and C-rich TP- 
AG B stars (bottom panel) , for the Rizzi et al. (in prep.) {solid line) 
and lDolphin et al.l boOSl. dotted line) SFHs. The grey histograms 
correspond to the stars actually observed. 



Classical models of stellar evolutior|f| have a strong difficulty 
in reproducing such low mass values for the progenitors of 
carbo n stars (see e.g. lHerwiell2005l : IStancliffe. Izzard fc Tout! 
l2005l ). 



7.4 Luminosity functions 

A simple comparison of Fig. [3] and [11] reveals that simulated 
C stars have about the same Ka magnitudes the observed 
ones. This is confirmed by the bottom panel of Fig |12l which 
shows the mean C star luminosity function (as derived from 
a total of 50 simulations considering the SFR errors) as com- 
pared to the data. A KS test indicates a 59 % probability 
that the observed distribution is drawn by the predicted one, 
in the case of Rizzi et al. (in prep.) SFR, a nd a 56 % prob- 
ability in the case of iDolphin et al.l (|2005l ). Also, we note 
that the presence of a fraction of C stars located below the 
TRGB is supported by present models. 

The top panel of Fig [T2] shows the comparison for O- 
rich giants, limited to the Ks = 15.75 interval so as to avoid 
RGB stars. Here, it is evident that the models predict a 
distribution extended to brighter magnitudes than observed. 
The KS probability that the observed distribution is drawn 
by the predicted one is very low, i.e. lo wer than 2 % for both 
Rizzi's et al. and lDolphin et al.1 (|2005h SFRs. 

The discrepancy of O-rich giants is statistically signifi- 
cant. The comparison with theoretical isochrones of Fig. [TT] 
suggests that the bright stars missing in the data can be 
identified either with young early-AGB stars belonging to 

^ For classical models we mean models following the 
Schwarzschild criterion for convective borders, and solving 
separately for the stellar structure and nucleosynthesis. 



populations younger than ^ 3 Gyr, or to the bright sec- 
tion of the TP-AGB for populations older than 7 Gyr. We 
have verified that eliminating any SFR younger than 4 Gyr 
from the models, the problem persists. Most of the excess of 
bright TP-AGB stars is produced by the old populations. 

Moreover, the discrepancy in the bright part of the 
LF is likely related to the excess of O-rich giants that we 
find in the simulations. In TP-AGB models without third 
dredge-up, the total lifetime depends essentially on the mass 
loss efficiency, and particularly on the critical region of 
the luminosity-temperature plane that triggers a superwind 
phase where most of the stellar envelope is lost. If the super- 
wind phase is delayed, both the lifetime and the luminosity 
excursion of the TP-AGB phase would be overestimated. 
Indeed, anticipating the superwind phase in O-rich models 
of low mass and metallicity may constitute an interesting 
solution to the discrepancy we find. 



8 SUMMARY AND CONCLUSIONS 

We have presented near-infrared JHKs photometry of a 
13!6 X 13'6 field centred on Leo II dSph, obtained with the 
new wide field imager WFCAM mounted at the UKIRT tele- 
scope. Our data cover most of the extension of Leo II dSph, 
and are complemented by optical data obtained with the 
EMMI camera at the ESO NTT telescope. 

The good statistics of our database, together with the 
wide colour baseline, allowed a precise determination of the 
distance and metallicity of Leo II. We derived a distance 
modulus (m - M)o = 21.68 ± 0.11 from the J, H, and Ks 
band magnitudes of the TRGB. This is in agreement with 
optical results, confirming the reliability of our NIR meth- 
ods. 

The V — Ks colours of RGB stars were used to de- 
rive the metallicity distribution of the stellar populations 
of Leo II. Using RGB fiducial lines of GCs as templates, we 
measured a mean metallicity [M/H] = —1.74. Since the bulk 
of the stellar population of Leo II is relatively old, we esti- 
mated that the population correction to be applied to this 
value is modest. Assuming a mean age of 9 Gyr yielded a 
correction of 0.10 dex, from which the age-corrected metal- 
licity is [M/H] = —1.64. Our measurement is in excellent 
agreement with recent spectroscopic results (jBosler et al.l 
l2007l : iKoch et al.|[2007l ) . A direct comparison between spec- 
troscopic and photometric metalli cities of ind i vidual stars 
suggests that the ages derived by iKoch et al.l (|2007l ) may 
be underestimated. Indeed, older mean stellar ages would 
be in better agreement with the SFHs obtained f rom HST 
photometry l|Hernandez et al.ll2000l : lDolphinll2o"o3 ). 

We also used our NIR data to define the properties of 
a nearly complete sample of AGB stars in Leo II dSph. By 
selecting AGB stars in the NIR two-colour diagram, we were 
able to discriminate the C-rich from O-rich stellar popula- 
tions. Foreground Milky Way stars are also easily separated 
from Leo II AGB members. Our NIR photometry was cross- 
identified with_^reviou^ of AGB stars, in particular 
C stars (lAzzopardi et al] 1985 ). One of the 7 carbon star 
listed bv lAzzopardi et al.l ( 1985 ) has anomalous colours, and 
visual inspection of our images confirms that it is a back- 
ground galaxy. No indication for additional C stars above 
the TRGB was found in our analysis, therefore we conclude 
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that the remaining 6 stars represent the complete population 
of C stars in Leo II within the area covered by our observa- 
tions. Using our colour selection, we provide the first sample 
of O-rich stars in Leo II above the TRGB, with negligible 
contamination from foreground Milky Way stars. 

Our Leo II observations were modeled via simulations 
based on the HST-derived SFHs, and using the most up- 
dated isochrones. The comparison between data and simu- 
lations has evidenced both successful and discrepant points, 
which are all potentially important for the calibration of 
AGB star models at low metallicity. 

With respect to the O-rich TP- AGB stars, the most im- 
portant discrepancy consists in a predicted over-estimation 
of their number and mean _K"3- band luminositi e s, as com- 
pared to the data. Interestingly, IWiUiams et al] l|2007l ) find 
sim ilar indications of a n excess of AGB stars in the old- 
est iGirardi et al] (|2000l ') isochrones while fitting RGB and 
AGB stars in the Virgo intracluster stars observed with 
HST/ACS. A possible solution to this problem could be an 
increase of mass loss efficiency in O-rich TP-AGB models 
of low metallicity and low mass. This will be explored in 
forthcoming work. Regarding this point, we note that: 

(i) The obscuration of AGB stars by circumstellar dust 
could be much more efficient than assumed here for low 
metallicity, and contribute to the solution of this problem. 
Indeed, recent Spitzer/IRAC obs ervations of the dwarf irreg- 
ular galaxies WLM and IC 1613 (|jackson et al.ll2007al lbh in- 
dicate a very high fraction of optically-obscured AGB stars, 
of about 40 %. In order to reduce the discrepancies in the 
LPs for Leo II, dust obscuration should be affecting mainly 
the O-rich AGB stars of higher luminosities. 

(ii) This discrepancy could still be reduced by adopting 
alternative SFHs. It would be again desirable to have im- 
proved derivations of the SFH in dwarf galaxies, based on 
HST observations co vering larger areas than those used by 
iDolphin et al.l (|2005h and Rizzi et al. (in prep.). Progress in 
this sense is expected as the result of o ngoing HST surveys 
and legacy programs on dwarf galaxies (|Dolphin et al.ll2005l : 
lDalcanton|[200d : lGallartll200'<i ). 

The above-mentioned discrepancy is still based on too a 
small sample of observed AGB stars. In forthcoming papers 
we will extend the comparison to other dwarf galaxies in the 
Local Group, in order to increase the statistical significance 
of our conclusions. 

With respect to the C-rich TP-AGB stars, the compar- 
ison between data and simulations is satisfactory, in terms 
of location in the colour-magnitude diagram, counts, and 
luminosity functions. We derive the following important in- 
dications: 

(i) At low metallicity the minimum mass for a star to 
become a C star can be as low as ~ 1 Mq , which sets an 
important constraint to the treatment of the third dredge-up 
in AGB stellar models. 

(ii) Low-mass C stars produced by the TP-AGB evolution 
of single stars (i.e. not belonging to binary systems), are 
expected to populate the O-rich sequence of giants below the 
TRGB when they are in the long-lived low-lumin o sity d ips 
driven by thermal pulses iBoothrovd fc SackmannI (|l988h . 

(iii) The higher effective temperatures of C stars with 
bluer J — A's colours can be a consequence of (a) less efficient 



molecular formation and opacity at lower metallicity, and/or 
(b) warming of their Haya shi line during the low-l uminosity 
stages of pulse cycles (see iMarigo fc Girardil [20071 '). 

From our analysis, we conclude that constraints to the 
TP-AGB models can be obtained from dwarf galaxies with 
known SFHs, provided that the numbers of AGB stars are 
significant and that the SFH is known with sufficient accu- 
racy. 
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